Contribution of weak localization to non local transport 
at normal metal / superconductor double interfaces 
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In connection with a recent experiment [Russo et al, Phys. Rev. Lett. 95, 027002 (2005)], 
we investigate the effect of weak localization on non local transport in normal metal / insulator / 
superconductor / insulator / normal metal (NISIN) trilayers, with extended interfaces. The negative 
weak localization contribution to the crossed resistance can exceed in absolute value the positive 
elastic cotunneling contribution if the normal metal phase coherence length or the energy are large 
enough. 

PACS numbers: 74.50.+r,74.78.Na,74.78.Fk 
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I. INTRODUCTION 

The manipulation of correlated pairs of electrons in 
solid state devices has aroused a considerable interest re- 
cently. The goals of this line of research are the realiza- 
tion of sources of entangled pairs of electrons for quantum 
information, and the realization of fundamental tests of 
quantum mechanics 1,2 . Correlated pairs of electrons can 
be extracted from a superconductor by Andreev reflec- 
tion, with extended or localized interfaces between super- 
conductors and normal metals or ferromagnets^^*SiL2iS, 
and in Josephson junctions involving a double bridge be- 
tween two superconductors^. 

Charge is transported by Andreev reflection at a nor- 
mal metal / superconductor (NS) interface: an electron 
coming from the normal side is reflected as a hole while 
a Cooper pair is transmitted in the superconductor—. In 
the superconductor, Andreev reflection is mediated by 
an evanescent state of linear dimension set by the super- 
conducting coherence length £. In structures in which 
a superconductor is multiply connected to normal metal 
electrodes separated by a distance of order the An- 

dreev reflected hole in the spin-(— a) band can be trans- 
mitted in an electrode different from the one in which 
the incoming spin-cr electron propagates. This "non lo- 
cal" transmission in the electron-hole channel is called 
"crossed Andreev reflection" . Non local transmission in 
the electron-electron channel corresponds to "elastic co- 
tunneling" by which an electron is transmitted from one 
electrode to another while spin is conserved^. 

A schematic three-terminal device is represented on 
Fig. as well as the voltages V a and Vf, applied on the 
normal electrodes "a" and "b" respectively. The volt- 
age Vs on the superconductor is chosen as the reference 
voltage (Vs — 0). Non local transport is characterized 
by a current I a (Vb) circulating in electrode "a" in re- 
sponse to a voltage Vj, on electrode "b". It is supposed 
in addition that V a = 0: electrode "a" is grounded, like 
in experiments^ (see Fig. QJ. The crossed conductance 
is defined by Q a ,b(Vb) = dI a /dVb(Vb). A crossed con- 
ductance dominated by elastic cotunneling (crossed An- 
dreev reflection) is negative (positive)^ because of the 
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FIG. 1: (Color online.) Schematic representation of the tight- 
binding model of NISIN trilayer. The insulating layers are not 
shown on the figure. The aspect ratio is not to the scale of 
the experiment in Ref. [3 for which d = 15, 50, 200 nm, L is of 
order of 5 /xrn, and d! ~ 50 nm. 



opposite charges of the outgoing particle in elastic co- 
tunneling and crossed Andreev reflection. Lowest order 
perturbation theory in the tunnel amplitudes leads to 
Ga,b(Vb) = because the crossed Andreev reflection and 
elastic cotunneling channels have in this case an exactly 
opposite contribution to the crossed conductance once 
the average over disordeni&ii or over the Fermi oscilla- 
tions in multichannel ballistic systems^^ is taken into 
account. 

Three unexpected experimental features for the 
crossed conductance in a normal metal / insulator / su- 
perconductor / insulator / normal metal (NISIN) trilayer 
have been reported recently by Russo et alA i) The 
crossed conductance does not average to zero with nor- 
mal metals, in contradiction to the abovementioned pre- 
diction of lowest order perturbation theory in the tunnel 
amplitudes^. The order of magnitude of the experimen- 
tally observed crossed signal is far from being compatible 
with lowest order perturbation theory, ii) A magnetic 
field applied parallel to the interfaces suppresses the non 
local signal, suggesting a phase coherent process, iii) 
The sign of the crossed resistance^ crosses over from 
positive (the sign of elastic cotunneling) to negative (the 
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sign of crossed Andreev reflection) as the bias voltage V& 
increases, and the crossed signal disappears if the bias 
voltage exceeds the Thouless energy in the superconduc- 
tor. 

We show here that weak localization with extended in- 
terfaces leads to a positive crossed conductance, the sign 
of which is opposite to the sign of the dominant elastic 
cotunneling channel for localized interfaces^. The weak 
localization contribution to the crossed conductance be- 
comes important at large bias voltages and for a large 
phase coherence length in the normal metals. 

The article is organized as follows. Necessary prelim- 
inaries are given in Sec. [H] The weak localization con- 
tribution to crossed transport is discussed in Sec. II I II 
Concluding remarks are given in Sec. IIVI 



II. PRELIMINARIES 

A. Hamiltonians 

The normal metal electrodes are described by the tight 
binding Hamiltonian 



B. Green's functions 



(1) 



where k is the wave-vector, a = j, [ the projection of the 
electron spin on the spin quantization axis, and where 
e(k) = —2to[cos (k x a Q )+cos (k y a Q )+cos (k z ao)] is the dis- 
persion relation of free electrons on a cubic lattice, with 
to the bulk hoping amplitude, ao the distance between 
two neighboring sites, and k Xl k y and k z the projections 
of the electron wave-vector on the x, y and z axis (see 

Fig. n) . 

The superconductor is described by the BCS Hamilto- 



Wbcs = £ ( k ) c iL c k,<T+ A Y (>.•'• k.. + c -k,l c k,1 

k,CT k 

(2) 

with A the superconducting gap. 

Diagonal disorder is included^, with an elastic mean 

free path le in the normal electrodes of the NISIN tri- 

f S) 

layer, and l e in the superconducting electrode. A fi- 
nite coherence length l^ 1 ' in the normal electrodes is ac- 
counted for by adding phenomenologically an imaginary 
part to the energy. 

At the extended interface a-a, the tunnel Hamiltonian 
connecting the electrodes "a" and "a" takes the form 



(3) 



where the sites a n on the normal metal side of the inter- 
face correspond to the sites a n on the superconducting 
side of the interface (see Fig.^l, and a similar expression 
holds for the tunnel Hamiltonian at the interface b-j3. 



The fully dressed advanced and retarded equilibrium 
Nambu Green's function G A ' R (w) at energy w is first de- 
termined by solving the Dyson equation 



G^M, (4) 



where ® denotes a summation over the spatial indices, 
and g A ' R (uj) are the advanced and retarded Green's func- 
tions of the disconnected system (i. e. with t — in the 
tunnel Hamiltonian given by Eq. J3Jl). The self-energy 
corresponds to the couplings in the tunnel Hamiltonian 
given by Eq. © . The Green's functions are 4x4 matrices 
in the spin ® Nambu representation. The four compo- 
nents correspond to a spin-up electron, a spin-down hole, 
a spin-down electron and a spin-up hole. Because of spin 
rotation invariance, some elements of the 4x4 Green's 
functions are redundant. We work here in a 2 x 2 block in 
the sector S z = h/2, encoding the superconducting corre- 
lations among a spin- up electron (Nambu label "1") and 
a spin-down hole (Nambu label "2"). 

Once the fully dressed advanced and retarded Green's 
functions have been obtained, the Keldysh Green's 
function G +: ~ is determined by the Dyson-Keldysh 
equation 



21^22^23,24.25 



I + G 1 



t) ® g + '~ ® (/ + tt ® G A ) , (5) 



where g + '~ is the Keldysh Green's function of the iso- 
lated electrodes, and where the energy dependence of the 
Green's functions is omitted. 



C. Transport formula 



The current through the interface a-a is given by 



i an , an Gt'- a M (6) 



ia n ,a n G^ an {uj) <T 3 | 



where the trace is a summation over the two components 
of the Nambu Green's function, 173 is one of the Pauli 
matrices, the diagonal elements of which are (1, —1), and 
the sum over n runs over all sites at the interface a-a. 
As shown on Fig. ^ the symbols a n and a n in Eq. © 
label two corresponding sites at the interface, in the nor- 
mal electrode "a" and in the superconductor respectively. 
The two spin orientations are taken into account in the 
prefactor of Eq. © . 

The local conductance of a NIN interface is equal to 
(e 2 /h)T per channel, where T = 2n 2 t 2 p 2 N is the dimcn- 
sionless transmission coefficient in the tunnel limit, with 
Pat the normal density of statesSi. 

The non perturbative transport formula for the lo- 
cal current at a localized NIN interface was obtained in 
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FIG. 2: (Color online.) A weak localization diagram in a normal metal for which ri and r2 are within the elastic mean free 
path l e (a), and in a superconductor for which n and T2 are within £. (c) represents the same diagram at a NISIN double 
interface, equivalent to Fig.|3J:. In the NISIN case on (c) ri and T2 are within £ and ri, V2 are within l\ ' from the interfaces. 
In the NIN case on (c), ri and r2 are within the insulator correlation length. "A" and "R" stand for advanced and retarded. 



Ref. l2ll and generalized in Ref. to a localized NIS 
interface, and in Ref. to non local transport at a dou- 
ble ferromagnet / superconductor interface. We deduce 
from these references the exact expressions of the local 
conductance GmN(V) of a single NIN interface, of the 
Andreev conductance Gms(V) of a single NIS interface, 
and of the crossed conductance G a ,b(Vb) of a NISIN tri- 
layer with extended interfaces: 

g N iN(V)=8n 2 ^t 4 (7) 

j,k,l 

Gms(V) = 16vr 2 i 4 (8) 
h 

j,k,l 

Ga,b(V b ) = -87T 2 ^t 4 (9) 

xE Tr {^,a/3G^ A (eH)^, bi <T 3 Gg iQs (eH)}, 

j,k,l 

valid to all orders in the hopping amplitude. The crossed 
conductance Ga,b(Vb) is per conduction channel through 
the junction of area L x L, with a superconducting layer 
of thickness d (see Fig. QJ. The summations in Eqs. JSJ 
and © run over all pairs of sites at the interfaces, and 
the overline denotes disorder averaging. The density of 
states p% ua . = (g^ a . - g* taj )/2m connects the sites a* 
and a,j in electrode "a", and a similar definition holds 
for Pb. b .. The site aj (pi) in electrode "a" ("b") is con- 
nected by the tunnel Hamiltonian to the site aj (/?;) in 
the superconductor (see Fig. [3] for the NISIN interface) . 
In the case of the local conductance of a NIN junction, 
aj, Pk, Pi and ai belong to an insulating layer that has 
been inserted in between the two normal metals. The 
tunnel amplitude t in Eq. J7J connects in this case the 
normal metal to the insulator, while t in Eqs. © and © 
connects the normal metal to the superconductor. 



The fully dressed advanced and retarded equilibrium 
Nambu Green's functions G„. p k (u>) and G^ a .(w) at en- 
ergy oj in Eqs. (J7J, JSJ and © are expanded by using the 
Dyson equation given by Eq. J5J . We deduce from the ex- 
act Keldysh transport formula given by Eqs. (Q, © and 
© that the diagrams are connected, and that the prop- 
agators for the density of states are directly connected 
to at least one tunneling vertex. The other extremity 
of the density of states propagators is connected cither 
to a tunneling vertex, or to a disorder scattering vertex. 
The density of states, represented by wavy lines on the 
diagrams, connects to one advanced and one retarded 
Green's function as in Eqs. Q, and ©• 

A crossed conductance dominated by elastic cotunncl- 
ing (crossed Andreev reflection) is negative (positive) — . 
This can be seen most simply by making the Nambu la- 
bels explicit in Eq. @ and taking into account the signs 
in the a 3 matrices, and the global sign. 

The crossed conductance is expanded perturbatively in 
T 2 , where the normal local transmission coefficient T has 
already been defined, and is also expanded in the number 
of non local processes crossing the superconductor since 

D. Weak localization in a superconductor 

Weak localization in a superconductor was already in- 
vestigated in Ref. |2(| in connection with the determina- 
tion of the number density of superconducting electrons. 
A weak localization diagram in a bulk normal metal is 
shown on Fig. In this case the two points ri and r2 
are within the elastic mean free path li N \ In a bulk su- 
perconductor (see Fig. |2>) , the two points ri and r 2 are 
within the coherence length £ since the disorder average 

9ap9ai3 °f t w0 advanced Green's functions between the 
two sites a and /? at ri and r 2 in a superconductor is 
limited by £, not by l e . Similar diagrams for the con- 
ductance are introduced at a NISIN double interface in 
Sec. IIIII In this case the points ri and r 2 are separated 
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FIG. 3: (Color online.) The lowest order diagrams contribut- 
ing to crossed transport. We show (a) the vanishing diagram 
of order T 2 , (b) the vanishing weak localization diagram of 
order T 4 , (c) the first non vanishing diagram of order T , 
and (d) the weak localization diagram of order T 6 . The wavy 
lines correspond to the insertion of the normal metal density 
of states, and of a 03 matrix. 



by the superconductor thickness, of order £ (see Fig- Et) • 
At a NIN interface, the two points ri and V2 on differ- 
ent interfaces can be separated by a distance equal to 
the thickness of the insulator, comparable to the decay 
length induced by the charge gap of the insulator. 



in Sec. EIU 

The diagram of order T 4 on Fig. [3: takes a finite value, 
and is evaluated explicitely by summing over all possible 
Nambu labels, and over the different possibilities of in- 
serting the density of states and the 173 matrices (repre- 
sented by the wavy lines on the diagrams). This leads to 
the crossed conductance 

GaAVt) = 4K 2 t 8 p%J2(&(eV b )g 2 % R (eV b ) (10) 



5aS' A ( e ^)ff^' R (eV6) 



A 2 



when 



,19 



1.2 / \ 1.2/ \ 

Sa,/j( w )ff/3,aM 



-VI 



-VI 



n(R) exp 



n(R) exp 



2R \ 2lu 2 - A 2 



CM 

2R 



A 2 



A 2 



with R the distance between the sites a and (3, and with 
n(R) = (fc^i?) 2 in the ballistic limit, and n(R) — kpRle 
in the diffusive limit. The resulting term of order T 4 in 
the crossed conductance is given by 



r {ec) 



04) 



rp4 



A 2 - (eV b y 



■ exp 



in agreement with the expansion to order T 4 of the 
crossed conductance obtained for highly transparent lo- 
calized interfaces 19 . A summation over the pairs of sites 
a and at the two interfaces was carried out, giving rise 
to the prefactor £/li S) in Eq. JTTJ- 



III. NON LOCAL TRANSPORT IN A NISIN 
TRILAYER 



C. Weak localization diagram of order T 



Crossed conductance to order T 



Contribution of local excursions 



Now we evaluate the lowest order diagrams appear- 
ing in the crossed conductance of a NISIN trilayer. The 
crossed conductance due to the diagram of order T 2 (see 
Fig. [3^) vanishes because the contributions of the elas- 
tic cotunneling and crossed Andreev reflection channels 
are exactly oppositei^iiSiiSii 9 .. This can be seen also by 
evaluating the summation over the Nambu labels in the 

diagram on Fig. Eland using fl^'jjfl^ = S&JsS&a) where a 
and f3 are two points in the superconductor at a distance 
of order £. 



B. Crossed conductance to order T 4 

The first weak localization diagram of order T 4 in- 
volving four Green's functions in the superconductor is 
shown on Fig. \3]p. The corresponding crossed conduc- 
tance vanishes, as for the diagram of order T 2 discussed 



We consider now the weak localization diagram of or- 
der T 6 on Fig. I2JI, merging the features of the two di- 
agrams of order T A on Figs. |2t> and c, and shown in 
more details on Fig. 0] The long range propagation in 

the normal electrodes involves the diffusons g^ 'p' 1 da 'p' 1 

and 5„ '|' 2 <?^'|' 2 in the particle-particle or hole-hole chan- 
nel (where a and (3 belong to the normal electrode), 

A 1 1 A 2 2 

as opposed to the diffuson g a 'p' g a 'g' in the particle- 
hole channel for local Andreev reflection at a single NIS 
interfac o 27 i 28 below the related Thouless energy. 

We first evaluate the part of the diagram on Fig. [3] in- 
volving local excursions at the NIS interfaces (see Fig.|SJ. 
Enumerating these diagrams shows that the two local 
excursions are attached to the same non local propaga- 
tion in the superconductor. We use the notation X^ A r , r , 

(A^:^, r ,) for the part of the diagram on Fig. [5] with two 
advanced (one advanced and one retarded) propagators 
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where & (V&) is given by Eq. (jTTJl and 



FIG. 4: (Color online.) Details of the insertion of the densi- 
ties of states and of the 0-3 matrices in the weak localization 
diagram on Fig. |3Ji. "A" and "R" stand for advanced and 
retarded. The diffusons are calculated in the ladder approxi- 
mation. The rungs of the ladders represented by the red lines 
represent schematically the diffuson in the ladder approxima- 
tion. 



in the superconductor, and the Nambu labels (r, r) and 
(t',t') at the extremities in the two normal electrodes. 
We find 



V A.A 
A ll,ll 



^11.22 — 



X 



A,R 
11,11 



V A.R _ 
^11,22 — 



2* 2 P % A 2 
H(R) A 2 - 



2 \2 



(12) 
(13) 



We recover Eq. (|ll|l for the non vanishing diagram of 
order T 4 on Fig.^t, proportional to — (X 1 i 11 — A^ 11 'f 2 ). 
The global minus sign in this expression can be found in 
Eq. (O, and the opposite signs for xf{^ and A^'f 2 are 
due to the matrices 03 in Eq. ijHJ . The diagrams with lo- 
cal excursions attached to the disorder average of two ad- 
vanced or two retarded Green's functions in the supercon- 
ductor lead to a vanishingly small crossed conductance, 



as it can be seen from the identity — (-Xix'11 — ^11 227 
deduced from Eq. (ft^)! . 



,A,A 



= 0, 



U' A',W 





(15) 



A 2 



(ev b y 



■ exp 



The sign of Q^" b l \eVb) given by Eq. (|l4fl is positive, as 
for crossed Andreev reflection. The factor (lp /le ) 4 is 
due to a a factor (l^ /li N ^) 2 associated to the diffusons 
in each of the quasi-two-dimensional normal layers [see 
Eq. i|A5(l ] in the limit q — 0. The constraint q = 
originates from the conservation of the component of the 
wave- vector parallel to the interface, due to the symmetry 
by translation parallel to the interfaces (see Appendix lB|l . 

It can be shown that the two diagrams of order T 6 in- 
volving a single diffuson in the superconductor and long 
range propagation in the normal metals are negligible 
because of the sum over the Nambu labels in one dia- 
gram, and because of the factor /li N ^) 2 in the other 

diagram, much smaller than (l^ /^e^') 4 for the weak 
localization diagram. 

The weak localization crossed conductance can be ex- 
panded systematically in powers of {l v /li N ^) 2 : 



L 



2n 



(16) 



The coefficients A n (T) are evaluated to leading order in 
T because of the small interface transparencies. Estimat- 
ing the higher order weak localization diagrams leads to 
A 3 ~ — T 8 , A 4 = 0, A 5 ~ T 10 , A 6 = 0,A 7 ~ -T 12 , ... 
The order of magnitude of the sum of the higher or- 
der contributions can be obtained from summing the 
corresponding geometric series for T{l v /li N) ) 2 > 1: 



E, 



A 



is of order T 2 



(l v /li N) )\ 



much 



smaller than A 2 {T) , therefore justifying why 

we based our discussion on the first two terms Aq(T) ~ 
T 4 and A 2 (T) ~ T 6 . 



E. 



Relation to experiments 



D. Weak localization crossed conductance 

The weak localization contribution to the crossed con- 
ductance corresponding to the diagram on Fig.^is given 
by 

G%HeV b ) = -G { a f(V b )F (j, ^, l -jpj , (14) 



1. Determination of the parameters 

The number of channels N c ^ for a contact with a three- 
dimensional metal is 7V c h = S/Xp, where S is the junc- 
tion area and Xf the Fermi wave-length. The normal 
layers have a dimension L x L x d', with L ~ 5 //m and a 
thickness d! ~ 50 nm (see Fig. QJ. The number of chan- 
nels in the quasi-two-dimensional geometry is obtained 
be neglecting disorder (the elastic mean free path in the 
experiment is limited by scattering on the normal film 
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FIG. 5: (Color online.) Representation of the part of the diagram of order T 6 (see Fig. [3J with local excursions. "1" and 
"2" correspond to the Nambu labels, (a), (b), (c) and (d) contribute respectively to Xj^'^, X^^ 2 , and X x {^ 2 . The 

diagrams having the topology of (e) and (f) lead to a vanishingly small crossed conductance because of the matrices 63 and 
the trace over the Nambu labels. 



boundaries), and by evaluating the area of the Fermi 
surface with discrete wave-vectors in the direction per- 
pendicular to the interface, leading to iV c h = Ld! j\ 2 F ~ 
6 x 10 5 . The normal transparency T can be obtained from 

the local conductance in the normal state 7Z 



(loc 
N 



K 



1 e 2 
~ 2N ch T— 

(loc) cn /j 

N 



2Q: 



(17) 



leading to T ~ 10~ 2 . These values T ~ 10" 2 and N ch ~ 
6 x 10 5 are compatible with the local Andreev resistance 
at zero bias of about 100 f2, being an upper bound to 
the zero temperature Andreev resistance TZ^ C ^ given by 
jy^jioc) = 2(e 2 /h)T 2 N ch . Russo et al& estimate T ~ 
10 -5 from N c h = S/\ 2 F for a three dimensional metal. 
The possible dependence of on dl that we consider 
here can be probed experimentally by determining how 
the crossed resistance depends on the thickness of the 
normal layers. 



2. Crossed resistance 

The crossed resistance matrix measured 
experimentally 9 is the inverse of the crossed con- 
ductance matrix. The off-diagonal element of the 
crossed resistance matrix is given by 



n^f(v b ) = 



Ga£\v b )g£™ } (v b ) - G y :T'(VbKT'(Vb) 

(18) 

where Q%?\V h ) = Q [ b f{V b ) = g£?\V b ) is the local 
Andreev conductance. 

The elastic cotunneling crossed resistance correspond- 



,(tot) , 



,(tot) , 



,(tot) 



ing solely to the contribution of Eq. I|llfl is thus of order 



K$(V b = 0) 



4A ch e 2 \i(S) 
(eV b ) 2 



A 2 



exp(-2d/£), (19) 



having an order of magnitude compatible with the 
experiment^. The elastic cotunneling crossed resistance 
is independent on T for tunnel interfaces. The elastic 
cotunneling crossed resistance T^ b {y b — 0) being in- 
versely proportional to A c h, with N c ^ proportional to d' 
(see Sec. IIII E Q , it is expected that the crossed resis- 
tance would decrease if the normal metal layer thickness 
d! increases, as anticipated in Sec. IIII E II 

Now, the weak localization crossed resistance is equal 

to 



= -n 



(ec) 



'd V b 



6 U' A 'il N) 



Jm ■ (20) 



The voltage dependence of the total crossed resistance 
K a , b (V b ) = ll { ™ b \V b ) + n^ b l) (V b ) is shown on Fig. © 

for different values of l^/li. We obtain a change 
of sign from a positive (elastic cotunneling dominated) 
to a negative (weak localization dominated) crossed re- 
sistance as the voltage increases (see Fig. |SJ. For suf- 
ficiently large values of the phase coherence length, we 

have F (d/£, 0, /li N ^ > 1 as discussed above, so that 

the crossed resistance is negative for all values of the bias 
voltage. 

The perturbative crossed resistance on Fig. tends to 
a finite value in the limit eV b = A. The determination 
of the crossed conductance around eV b — A is examined 
in Ref. [19| for localized interfaces with arbitrary inter- 
face transparencies. The non perturbative crossed con- 
ductance tends to zero for eV b = A, but the cross-over 
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FIG. 6: (Color online.) Voltage dependence of the crossed 
resistance (in Ohms), given by Eq. I|180 . in which the crossed 
conductance is the sum of the elastic cotunneling [see Eq. 
and weak localization [see Eq. 1141 1 contributions. We used 



the parameters d = 20 nm, £ = 15 nm, I 



(N) 



3 nm, T = HT 2 
top to bottom: 



JVch = 6 x 10 5 



50 nm, I 

,(JV) 



0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1 ^m. The 



points are theoretical, and the lines are a guide to the eye. 



occurs within an energy window that becomes very small 
for tunnel interfaces. A crossed current related to out- 
of-equilibrium quasiparticle populations (not described 
here) is expected for eVb > A. 

As it is visible on Fig.[fjl the characteristic voltage scale 
in the bias voltage dependence of the crossed resistance 
is the superconducting gap, not the normal state super- 
conductor Thouless energy obtained in experiments. At 
the present stage, we do not find a plausible explanation 
of this experimental observation. 



IV. CONCLUSIONS 

We have calculated the weak localization contribution 
to non local transport in NISIN trilayers with extended 
interfaces and a sufficient phase coherence length in the 
normal electrodes. We find a change of sign in the crossed 
resistance between elastic cotunneling at low voltages and 
weak localization at higher voltages. The weak local- 
ization contribution can dominate for all voltages if the 
phase coherence length is large enough. The weak local- 
ization crossed conductance dephases in an applied mag- 
netic field, but not the elastic cotunneling contribution. 
The appearance of a voltage scale related to the super- 
conductor Thouless energy is left as an important open 
question. 

The author thanks H. Courtois, A. Morpurgo, and S. 
Russo for many useful suggestions and comments. Spe- 
cial thanks to D. Feinberg and F. Pistolesi for numerous 
fruitful discussions. 



APPENDIX A: DIFFUSION PROBABILITIES IN 
A NORMAL METAL 

1. Diffusion in the electron-electron channel 

We first discuss briefly the normal metal diffuson in the 
ladder approximation. In the Born approximation, the 
elastic scattering time T\ t \ in the "11" Green's function 
given by 



9i i (k,w) 



(Al) 



is defined by kpri^v 2 /{^Trep) = 1, with rx,i = r [l + 
w/(2ejr)], where tq is the elastic scattering time at the 
Fermi level, v is the amplitude of the microscopic im- 
purity scattering potential, fop is the Fermi wave-vector, 
and 6f the Fermi energy. 

Using contour integration, and the identity 



F. Magnetic field dependence 

The experimental crossed signal is suppressed by a 
magnetic field parallel to the layers^. The theoretical 
weak localization crossed signal is also suppressed by a 
magnetic field because the corresponding diagram de- 
phases in an applied magnetic field. The cross-over mag- 
netic field B,f for the suppression of the weak localization 
crossed conductance corresponds to one superconducting 
flux quantum </>o enclosed in the area of the diagram, 
compatible with experiments 9 . However, in experiments, 
the crossed resistance is suppressed by a magnetic field 
in the entire voltage range. The present model does not 
explain the dephasing of the elastic cotunneling contri- 
bution. 



(A2) 

with k = ±|k| and /(k) a function of k, we obtain easily 
/ (^^(^(k + q) (A3) 

where D is the diffusion constant. The Fourier transform 
P N (q) of the diffusion probability V N (r), given by 

^-j-^Gf^k^G^ (k + q, W ), (A4) 
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where the overline is a disorder averaging and to is the 
bulk hopping integral (see Eq. QJ), is obtained by sum- 
ming the ladder series, leading to 



metals are quasi-two-dimensional so that q± — 0, where 
q± is the projection of q on the normal to the interfaces. 



?> iV (q) 



Dt , 



(A5) 



In this expression V N (q) has no dimension. Its Fourier 
transform V (R) is such that V N (R)cZR has no dimen- 
sion, as expected for a probability. 



APPENDIX B: DOUBLE NS INTERFACE 

We provide now explanations to the factor (Z^, /Ze ) 4 
appearing in Eq. l|15fl for quasi-two-dimensional normal 
electrodes. Because of the conservation of the compo- 
nent of the wave- vector parallel to the interface at each 
tunneling vertex, and because of the form of the dia- 
gram, the diagram on Fig. 0] is evaluated at qn = 0, 
where q|| is the projection of q (see Appendix |A|) on a 
plane parallel to the interfaces. By momentum conser- 
vation, a finite value of q|| is transformed in — q|| after 
traversing the entire diagram, so that = — = 0. 
The resulting crossed conductance is thus proportional 
to [V N {q ]l = 0)] , where the square is due to the corre- 
lated diffusion in the two normal electrodes, and where 
the diffusion probability T' N (q\\) is given by Eq. O. 

This lead to the factor ~ (/£ ) 4 since the normal 



APPENDIX C: BALLISTIC NISIN TRILAYER 
WITH ATOMIC THICKNESS 



We consider in this Appendix a ballistic NISIN trilayer 
with atomic thickness^^S in which the three electrodes 
are two-dimensional, and find a factor (l^ /li N ^) 3 , sim- 
ilar to the factor (Z^/Ze ) 4 factor discussed previously 
in the diffusive limit. 

The transport formula given by Eq. @ is Fourier 
transformed, to obtain 



Tr 



d 2 k 

N 



(CI) 



(k,e!4)}, 



and the Dyson equations for G A ' R (k) are also Fourier 
transformed. The normal metal ballistic Green's func- 
tion is given by g^(k, lu) = 1/[lu — £ + irj\, where £ is the 
kinetic energy with respect to the Fermi level, and where 



hvpl 



(AO 



the broadening parameter r\ is given by r\ 
Eq. i|Cl|) is then expanded diagrammatically. A dia- 
gram similar to the one on Fig. 0] leads to a factor 
r?- 3 - (hv F l ( v N) f. 
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